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Superconducting thin-film rings smaller than the film penetration depth (the Pearl length) are 
considered. The current distribution, magnetic moment, and thermodynamic potential T(H,N,v) 
for a flat, washer-shaped annular ring in a uniform applied field H perpendicular to the film are 
solved analytically within the London approach for a state with winding number N and a vortex at 
radius v between the inner and outer radii. 

PACS numbers: 74.78.-w, 74.78. Na 



I. INTRODUCTION 



Small flat rings made of thin superconducting films 
are of interest for a variety of mesoscopic experiments^* 2 * 
SQUID-type devices is an example; another one is the 
study of fluxoid dynamics in rings, as discussed by Kirtlcy 
et al.i The basic physics of such rings is governed by flux 
quantization, as in the Little-Parks experiment; see, e.g., 
Ref. or later publication a 5 : 6 : 7 : 8 based on the Ginzburg- 
Landau theory. 



There is a resurgence of experimental interest in this 
subject: The interacting dipole moments in an array of 
superconducting rings provide a model system for mag- 
netism in Ising antiferromagnetsJ* 2 " Moreover, there is 
considerable interest in quantum coherence effects in su- 
perconducting rings and their arrays for potential appli- 
cations in quantum computingAifiiAi 



Many quantitative details specific to the thin-film ge- 
ometry can be treated within the London approach, 
which is not bound by the strict temperature restriction 
of Ginzburg-Landau models. In the thin-film limit, for 
which the London penetration depth A obeys A>d, the 
film thickness, the fields and currents are governed by 
the Pearl length A = 2\ 2 /dH As we show below, when 
the inner and outer radii a and b of the annular ring 
are smaller than A, it is possible to obtain analytic so- 
lutions for the energy of the ring in a uniform magnetic 
field H applied perpendicular to the ring plane with a 
vortex sitting in an arbitrary position at the annular re- 
gion between a and b. The motion of vortices between 
a and b provides the means for the ring to switch be- 
tween discrete states with different winding numbers N 
defined below. Study of these transitions is relevant for 
understanding the telegraph noise observed in multiply 
connected mesoscopic superconducting devices in general 
and in thin-film rings in particular^ 



II. CURRENT DISTRIBUTION 

Let us consider a small thin-film ring of thickness d <C 
A situated in the plane z — with inner and outer radii 
a and b, where b is much smaller than the Pearl length 
A = 2A 2 jd. The ring is in a uniform applied field H 
perpendicular to the film plane; it may contain a vortex 
(or antivortex), the position of which can be taken as 
x = v, y = (with the origin at the ring center and 
a < v < b). The London equations for the local magnetic 
field b in the film interior readi*^ 

47rA 2 

b H curlj = ±4>oz5(x — v,y). (1) 



Here, j is the current density, z is the unit vector perpen- 
dicular to the film plane, and cj>o = irhc/\e\ is the value 
of the flux quantum. The upper sign holds for a vortex 
whereas the lower one is for an antivortex, the convention 
retained throughout this paper. Averaging this over the 
thickness d, one obtains 
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cm\ z g = ±(j>ob~(r - v), 



(2) 



where g(r) is the sheet current density, r = (x,y), and 
v = (i>,0). Equations QJ and are valid everywhere 
in the film except within a distance of the order of the 
coherence length £ (vortex core) from the vortex or an- 
tivortex axis, where the London equation Q no longer 
holds. 

The distribution g(r) can be found by solving Eq. J5J), 
combined with the continuity equation and the Biot- 
Savart integral that relates the field b z to the surface 
current: 



divg = 0, b z (r) 



where R = r — r' 



J [g(r') x R/cR 3 } z d 2 r' + H, 



(3) 
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A. The stream function 

In principle, Eqs. J2J and determine the current 
distribution. To solve these equations for the general case 
is a difficult task, even for a disk4^ However, for small 
samples, as in our case for which b <C A, the problem can 
be solved^ To this end, let us introduce a scalar stream 
function G(r), such that 



g = curl(G£). 



(4) 



The first of Eqs. J3J) is then satisfied. It is easily seen that 
the contours G(x, y) = const coincide with the current 
streamlines. Substituting Eq. Q into (J2J, we obtain: 
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lV 2 G = T0o<5(r- - v) + b z 



(5) 



The radial component of the current at the ring edges 
must be zero; this means that the values of G(x,y) at 
the edges are constants: 

G(r = a,cp)=G a , G(r = b,<p) = G b , (6) 

where (r, tp) are polar coordinates. Since g v = —dG/dr, 
the total counterclockwise current (in the p> direction) 
around the ring is 



I — G a — Gb 



(7) 



The self-field of the ring currents within the ring can be 
estimated using the Biot-Savart law: J cPrg/cR 2 ~ g/c. 
Substitution of b z ~ H + g/c into Eq. (J2J reveals that the 
self- field can be disregarded because i? C A; i.e., we can 
set b z = H in Eq. (J5J: 
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ly 2 G = Tfo${r - v) + H. 



(8) 



Since this equation is linear, we can look for a solution 
of the form G = G v + Gh , such that G v satisfies 



C0 O 



and 



(9) 



(10) 



One can say that G v describes the currents due to the 
vortex or antivortex, whereas Gh is the response to the 
applied field. 

The boundary conditions (|SJ| are imposed on the sum 
G v + Gh ■ It is convenient to require that 



and 



G v (a) = G v (b) = 0, 



Gjj(a) — G a , G H (b) = G b 



(11) 



(12) 



In a uniform field H, Gh can be taken as cylindri- 
cally symmetric. Aside from an unimportant additive 
constant, 



G H (r) = 



— r 2 +G ln- 
87rA a 



(13) 



where the constant Go is expressed in terms of the total 
current J7J): 



cH(b 2 - a 2 ) b 
1 8^V G ° ln a 



(14) 



To evaluate I we use the London equation in the form 

3 = -^r(w + ^). (15) 
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Here 9 is the order parameter phase, the topology of 
which plays a major role in our problem. 

When the ring is traversed around a circle of radius 
r in the positive direction of the azimuth ip, the phase 
changes by —2irN, where N is an integer which is com- 
monly called the winding number or the vorticity. If there 
are no vortices in the annulus, the integer N is the same 
for any contour within the annulus and we consider the 
state as having winding number N . In other words, the 
state of the system is characterized by the integer N and 
the continuous variable H . In zero field, N > cor- 
responds to positive currents g v and positive magnetic 
moments fi z . 

In the presence of a vortex, however, the situation is 
different. For contours encircling the ring's hole, the 
winding number is N at contours that do not include 
the vortex position; the number is N + 1 for those that 
do. The state is now characterized by the vortex position 
v in addition to N and H. In this paper, when the state 
of the ring with a vortex is characterized by the variables 
N,H, and v, it is implied that the integer N describes 
the phase topology on contours that do not include the 
vortex. The generalization to antivortices is obvious. 

Coming back to evaluation of the total current, we 
integrate I — J a g<p(T, tp)dr over tp to get: 

I=f g v (r)dr, (16) 

J a 

where the azimuthal average of g v (r, p) is 



2tt 



for r < v, and 



dp afro 



5 ^ 
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N ttH . 
l~ r . (I 7 ) 



c<fr Q (N±l irH \ 

We now readily evaluate the total current and the con- 
stant Gq: 



Go 



ccfro 



4tt 2 A 



\n(b/v) 



ln(6/a) 



(19) 
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Note that Go, which determines the part Gh of the 
stream function, depends on whether or not a vortex is 
present at the ring and on its position. In other words, 
the current distribution is not a simple superposition of 
the currents generated by a vortex in zero field and of 
those existing in the state (N, H) in the vortex absence. 



B. Electrostatic analogy and exact solution 

To find the solution of Eq. for the vortex-generated 
stream function G v subject to the boundary conditions 
(|II|I . we observe that the problem is equivalent to the 
two-dimensional one for the electrostatic potential gen- 
erated by a line charge ±c0o/87r 2 A at the point v sit- 
uated between two coaxial grounded metallic cylinders 
with radii a and b. The necessary conformal mapping 
procedure is given in Ref.E 



G v (w,v) 
A{w, v) 
B(v) 



4"7T 2 A 



In 



A(w, v) 



B(v) 



w 



re 



up 



cn[27 ln(v/a), m] cn[7ln(«/iy), m] 
sn[27 ln(u/a), m] 
cn[27 ln(i;/a), to] 



sn[7 ]u(v/w), to] 
cn[7 hx(v/a), m] 
i[2"/ hx(v/a), to] sn[7ln(u/a), m] 



(20) 



where cn(a;),sn(x) are the Jacobi elliptic functions, 7 = 
K(m)/ln(6/a) with K(m) being the complete elliptic in- 
tegral (in the notation of Ref . \l% . and the parameter m 
is chosen to satisfy 



K(l - to) ln(6/a) = ttK(to). 



(21) 



As an example, we find that for b/a = 2 (such rings were 
studied in Ref. 3), m = 1.048 x 10~ 5 , and 7 = 2.266. 

For rings with 1 < b/a < 2, we have < \n(b/a) < 1. 
Solving Eq. I|21|) numerically, one can see that to <C 1. 
Expanding functions K(to) and K(l — to) for small to, 
one obtains with a high accuracy: 



m = 16 exp 



ln(6/a) 



21n(6/a) 



One can now set to = in Eq. 1)20(1 to obtain 



G v 



ccpo J sin[7rln(vw/a 2 )/21n(6/a)] 
i - — tt R c \ .1: • 
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[it ln(w/w)/ ln(6/a)] 



(22) 



(23) 



At the vortex position r — v, ip — 0, this function 
diverges logarithmically. One can find G v (v) [and the 
self-energy e v = 0o|G,j(i>)|/2c] by introducing the stan- 
dard cutoff at a distance £ from the vortex axis at v: 



8tt 2 A 



In 



2vhi(b/a) . 7rln(u/a) 



ln(6/a) 



(24) 



The relative difference between this expression and the 
exact energy is less than 2xl0 -12 for6/a < 2>i& Formally, 




FIG. 1: Current streamlines in a ring with b/a = 2, N = 2, 
and vortex or antivortex position v/a = 1.2 in the normalized 
field h — H/Ho — 2.3 [see Eq. 139H 1. The streamlines are 
plotted as contours of G — G v + Gh = const. The left panel 
is for a vortex; the streamlines shown that encircle the hole 
have clockwise flow, whereas the flow around the vortex is 
counterclockwise. The right panel shows an antivortex; all 
streamlines shown have clockwise flow. 



the logarithmic factor in Eq. I|24|) goes to —00 if v — > a 
or v — y b. However, this expression fails when the vortex 
is within roughly £ of the inner or outer radius. In fact, 
e v of Eq. (|24|) becomes equal to zero at v a = a(l + £/2a) 
and Vb = b(l — £/2a). 

Thus, the problem of the current distribution in the 
ring is solved: the vortex-generated part G v of the stream 
function is given in Eqs. (|20|l . while Gh is determined by 
Eqs. 1(13(1 and 1(19(1 . The current streamlines are given by 
the contours G(r, ip) = G v (r, ip) + Gh(t) — const; Fig. 
shows examples of current streamlines for a vortex and 
an antivortex. 

It is worth mentioning that the same method based 
on application of conformal mapping to problems of the 
two-dimensional electrostatics can be utilized to obtain 
current distributions for n vortices equally spaced along 
a circle of radius a < v < b of the ringii& 



III. ENERGY 

A. Zero applied field 

Let us first consider the energy E(N, v) of the ring with 
a vortex at v in the state N in zero applied field. This 
energy is a sum of the magnetic and kinetic contributions. 
The magnetic part is^ 



1 

2c 
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■gd 2 r, (25) 



where we use Eq. (| 1 51) to express the vector potential A at 
the ring in terms of g and the phase 9. The supercurrent 
kinetic energy is the integral over the film volume of 
the quantity 27rA|j 2 /c 2 = it kg 2 / c 2 We find readily 



4 



that this energy is equal in value and opposite in sign to 
the term containing g 2 in Eq. I|25|) . so that 



E(N, v) = E m + E k = ~$!L J {W6 ■ g)d 2 r. 



(26) 



The integrand here can be further transformed in terms 
of the stream function G: 

V6> • cur\(Gz) = div(Gz x V0) + Gz ■ curlV0. (27) 

Substituting this into Eq. I|26|) . we use Gauss's theorem 
to evaluate the contribution of the first term: 



J div(Gz x V9)d 2 r 



-G{b,<p)— +G(a,<p) — 



= 2n[(N±l)G b -NG a }. 



dip (28) 



Here we have used the boundary conditions © and the 
phase change is —2ttN upon circling the inner radius and 
—2ir(N ± 1) the outer radius, provided a single vortex 
(antivortex) is present in the annulus. 

Integrating the second term of Eq. I|27|) . we use the 
basic quantization property of the phase [curl 2 V# = 
^2n5(r — v)], take the curl of Eq. l|15|l . and compare the 
result with Eq. J5J. Then, we obtain the magnetic and 
kinetic energy of the persistent currents in the ring: 



E(N,v) = ^-[G(v)-{N±l)G b + NG a \. (29) 
2c 

Note that G(v) = G v (v) + Gh(v), and the vortex self- 
energy is e v — 4> \G v {v)\/2-k (see Ref . [TEI) . Utilizing Eqs. 
(TT3)) and we obtain: 



E{N,v)=e v {v) + e 



N± 



ln(b/v) 



ln(6/a) 



where we introduce the energy scale 

4>lHb/a) 



eo 



8tt 2 A 



(30) 



(31) 



If v = b, i.e., if there is no vortex in the annulus, we 
have 



Eo = 7r(G a - G b )N 
2c 



2c 



-NI. 



(32) 



B. Magnetic moment 

By definition, the z directed magnetic moment readsi^ 



1 f 2n 7T 

V = 7T / drr 2 dipg v (r,(p) =- r^g^dr. 
2c J a Jo c J a 



(33) 



This is easily evaluated using Eqs. JTJJ, (fT5|) : 



» = +°<f-*)( N± * 



8ttA 



_^) _ A (34 ) 



It is worth noting that the magnetic moment is evaluated 
here without need to employ the explicit current distri- 
bution which could have been obtained from the known 
stream function G v + Gh- 

Setting v = b, we obtain the moment fio in the "pure" 
state N with no vortex at the ring: 



Mb 2 



8ttA 



-N ■ 



H 
16A 



(& 4 



(35) 



The result can be verified by direct calculation of fi using 
the current distribution in the absence of vortices. If we 
set v = a in (|34|l . we obtain no(N ± 1, H). 

It is instructive to observe that the exact moment H34|) 
in the presence of a vortex at v can be considered as a 
sum of the magnetic moment of a vortex-free ring with 
inner and outer radii a and v having winding number N 
and another ring with inner and outer radii v and b with 
vorticity N + 1. 



C. Free energy in field 

The thermodynamic potential T which is minimum 
in equilibrium at a given applied field is defined by a 
differential relation 



5T = n- SH = -liSH, 



(36) 



where fi(N, v, H) is the system magnetic moment^ Us- 
ing Eq. I|34|) , we readily obtain T by integrating Eq. I|36|) 
over H from to H: 



T ± =E(N,v) 



^H{b 2 -a 2 ) ( N± b 2 -v 2 



8nA 

+ ii^(6 4 -a 4 ), 
32A V h 



b 2 - a 2 



(37) 



where the zero-field energy E(N,v) is given in Eq. I|3U|I . 
After simple algebra we arrive at our main result, 



T±(N,v,H) = e„(u)+ e 



N± 



ln{b/v) 
\n(b/a) 



b 2 - v 2 i 



, (38) 



where we have introduced a dimensionless field h 
H/Hq, where 



2^0 ln(6/o) 



ir(b 2 - a 2 ) ' 

and a geometric factor \ given by 

_ b 2 /a 2 + 1 b 
X ~ b 2 /a 2 -l n a' 



(39) 



(40) 
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For a narrow ring b/a — l = ry^l, % = 1 + 0(i] 2 ). With 
increasing 6/a, x grows slowly: for 6/a = 2, x ~ 1-155, 
and x ~ ln(6/a) for large values of b/a. 

D. Vortex- free state 

For the vortex-free state, Eq. (|38fl yields: 

T (N,h) = e (^N 2 -2Nh + h 2 x). (41) 

It is seen that T$(Ni,h) = J 7 (N 2 ,h) at the field ft- = 
(Ni+N 2 )/2 for any ATi and N 2 . In particular, for ATj = N 
and i\T 2 = N + 1, the energies F (N, h) = T (N + l,h) 
at h = N + 1/2. In other words, at this field the system 
might be in either of the states N or N + 1 having the 
same energy. 

One can readily check that the thermodynamic poten- 
tial !Fo(N,h) of the vortex- free state coincides with the 
kinetic energy of the supercurrents in the ring. 

A transition from a state with vorticity N to one with 
vorticity N + 1 can happen when a vortex, carrying unit 
vorticity, enters at the outer radius, crosses the ring, and 
annihilates at the inner radius. Alternatively, such a 
transition can be accomplished, starting with initial vor- 
ticity N, when an antivortex, carrying vorticity -1, enters 
the annulus at the inner radius a, crosses the ring, and 
annihilates at the outer radius b, leaving behind vorticity 
N + 1 trapped in the ring. 

The field dependence of T§ for < A" < 4 is shown 
in Fig. [21 which illustrates the above features. It shows 
that in fields AT — 1/2 < h < N + 1/2, the minimum 
energy belongs to the state N. However, the first excited 
state is Af — 1 for N — 1/2 < h < N, whereas it is N + 1 
for N < h < N + 1/2. It should be also noted that in 
a given field h, the ground state winding number is the 
integer nearest to h. 

Although the ground state vorticity A~ changes with 
field, the energy difference between the lowest and the 
"first excited state" is, in fact, periodic in h with the 
period Ah =1. It is easy to check that the difference 
F (N, h)-F (N-l,h) = F (N+l, h+l)-T {N,h + l). 
This fact has implications for the transition probabilities 
from the ground state to the nearest excited state. 

It is worth observing that the energy (|38|l differs from 
that of the vortex- free state, Tq of Eq. I|41|) . by the vor- 
tex self-energy e v (v) and by the addition to N of the 
w-dependent terms that vary from unity to zero as the 
vortex moves from a to b. 



E. Potential barriers for vortices crossing the ring 

We are now in a position to evaluate the energy barrier 
that must be overcome for a vortex or an antivortex to 
move between the inner and outer radii of the ring. Since 
G v = at both a and b, we have at v = b: 

T+(N,b)=T-(N,b)=T (N), (42) 



3 
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FIG. 2: The energy JT of vortex-free states of a ring with 
b/a — 2 for vorticities TV from to 4 versus applied field. The 
energy and field units are defined in Eqs. 131H and 139H . Note 
that e.g. for 0.5 < h < 1.5, N = 1 corresponds to the ground 
state. The nearest "excited" state is N = for 0.5 < h < 1, 
whereas for 1 < h < 1.5 it is TV = 2. 

where the energy J-'o(N) of the vortex- free states is given 
by Eq. 1411) : the subscript "+" is for a vortex, whereas 
"-" stands for an antivortex (the field argument H is 
suppressed for brevity). In other words, the addition of 
either a vortex or an antivortex at the outer radius does 
not change the vorticity of the system. On the other 
hand, at v = a we have: 

T + (N,a)=T (N + l), T_(N, a) = T^(N - 1); (43) 

i.e., moving a vortex (antivortex) from b to a changes 
the vorticity by +1 (—1), an obvious consequence of the 
system topology. 

Starting from a vortex-free state N, the motion of a 
vortex between a and b can change the system energy in 
two ways. If a vortex moves in, starting from the outer 
radius b, the spatial dependence of the potential energy 
barrier against vortex entry is given by the difference 

r+(N,v)-F {N)=V£(N,v), (44) 

so that V+{N, b) = 0, while V+(N, a) = e (2N + l - 2h), 
which corresponds to the energy of transition between 
"pure" states N and AT + 1. On the other hand, if a 
vortex moves out, starting from the inner radius a, the 
spatial dependence of the potential energy barrier against 
vortex exit is given by the function 

V+ t (N,v)=F + (N-l,v)-F (N), (45) 

such that V+ t (N,a) = 0, while V+ ut {N,b) = e (2h - 
2A^ + 1), which corresponds to the energy of transition 
from the state Af to the state N — 1. 

When h = N + l/2, we obtain V+(N, a) = V+(N, b) = 
0, V+ t (N + l,b) = V+ t (N+l,a) = for any N. More- 
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over, the potential barriers V^(N, v) and V^ t (N + l,v) 
are identical. 

Similarly, the motion of an antivortex between a and 
b can change the system energy in two ways. If an an- 
tivortex moves in, starting from the outer radius b, the 
spatial dependence of the potential energy barrier against 
antivortex entry is given by the function 



VZ{N,v)=F-(N,v)-? {N), 



(46) 



so that Vr(N,b) = 0, while Vr(N, a) = e (2h- 2N + 1), 
which corresponds to the energy of transition from the 
state iV to the state N — 1. On the other hand, if an 
antivortex moves out, starting from the inner radius a, 
the spatial dependence of the potential energy barrier 
against vortex exit is given by the function 



V- t (N,v)=F-(N + l,v)-F (N), 



(47) 



such that V a ~ t (N, a) = 0, while V ~ t (N, b) = e (2N + 1 - 
2h), which corresponds to the energy of transition from 
the state N to the state N + 1. 

Note that V~ ut (N, v) differs from Vr (N + l,v) only by 
a constant: V^ t (N,v) - Vr(N + l,v) = F (N + 1) - 
T (N) = e (2N + l-2h). When this constant is zero, 
i.e. when h = N + 1/2, the barriers V~ ut (N — 1, v) and 
(N, v) become identical. 

We conclude the discussion of barriers by pointing 
out that although the total energies J-±(N,H,v) are 
quadratic in N and H, the barrier functions are linear 
in these variables. For example, we have for a vortex 
entry at b: 



V£(N,H,v)=e v (v) 
ln(b/v)~ 



(48) 



2N- 



ln(6/a) 



\n(b/v) 
ln(6/a) 



2h 



b 2 



IV. NARROW RINGS 

The energy (I38|l simplifies for narrow rings of width 
W = b — a <C a, i.e., for 

rj = W/a«l. (49) 

We obtain in the linear approximation in 77: 

T ± = e v (z) + e (N±z-h) 2 , (50) 

'° lnf^-sinTrz) , z=\ -. (51) 



4>t , /2a?7 . \ b-v 

—, 5— m — — sin ttz J , z = 

87r 2 A V Trf / b - a 



Here, < z < 1; the energy and field scales are 



87T^A 7TCT 



In particular, we have for the vortex- free state (2 = 0): 



fioV / at m2 _ ^o 7 ? Z'- 



Btt 2 A 



(AT - /i) 2 = 



8tt 2 A 



N - 



(53) 




FIG. 3: The energy barriers V^(N, H, v) for the vortex entry 
at the outer radius b calculated numerically according to Eq. 
1441 . The energy is given in units of (f) 2 /8n 2 A. The vortex 
position v is normalized on the inner radius a. The left panel: 
b/a = 1.2; the parameters from top to bottom are h = N = 
3, 110, 295, and 450; h is normalized to H of Eq. Pty . The 
right panel: b/a = 2; the parameters from top to bottom are 
h — N = 3, 10, 28, and 50. Note: the normalization field Ho 
for this set differs from the first. 



an expression similar to that for thin cylindersASS 

The potential barriers defined in Eqs. (|44fl - l|47[) are 
now easily evaluated. We focus here on for a vortex, 
which reads in linear approximation in 77: 



F 



+ e z (2N - 2h + z) . 



(54) 



One readily verifies that this potential reaches maximum 
at z m = 1/2 + 0(77), i.e., at v m = (b + o)/2. 

To study the behavior of the potential V^(N,H,v) 
near its maximum in increasing fields, one should go to 
higher order terms in the small 77. To avoid cumber- 
some algebra, we use the fact that in the ground state 
\h — N\ < 1/2 so that in large fields we can set A^ « h 
and study the behavior of the function Vjl(h,h, v). Of 
a particular interest is the curvature of the potential 
barrier at v = v m . Numerical experimentation shows 
that in large fields, the maximum of is situated close 
to v m = (b + a)/2. Differentiating twice the function 
V^(h,h,v) with respect to v (this is easily done using 
Mathematica) and setting v = v m = a + W/2, we find 
that the curvature turns zero at 



H « H-, = 



71" 00 

4IF 2 



(55) 



For H > Hi, the potential barrier V^{h, h, v m ) acquires 
a local minimum. Examples are shown in Fig. [21 A 
vortex at this minimum is in a metastable state provided 
V+(h,h,v m )>0. 

With increasing field, the depth of the minimum in- 
creases and at a field that may be called the low critical 
field, H c i, V^(h,h,v m ) becomes equal to zero. This is 
the minimum field at which vortices can nucleate at the 
ring and stay there in stable equilibrium. We can esti- 
mate this field by setting V^{h,h,v m ) = in Eq. lj4l?)) . 
For narrow rings we obtain: 



20o . 2W 
In — 



irW 2 



(56) 
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Clearly, with the further field increase, the barriers 
on both sides of the minimum are suppressed, while the 
points v where V£ = are pushed towards the annulus 
edges at a and b. The critical field H* at which the "edge 
barrier" near b disappears can be estimated by setting 
Vfo(h,h,vo) = in Eq. I|48(l and considering the limit 
vq w b — £. After straightforward algebra we obtain: 



H* 



6 In 2 



ln(6/o) 



tt£& 21n(6/a) - 1 + a 2 /b 2 

For narrow rings, this reduces to 

_ 0o In 2 
~ 2tt£W ' 

whereas for rings with a small hole, a -C 6, 
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1 



21n(6/a) 



(57) 



(58) 



(59) 



It is of interest to note that the field H% of Eq. I|55|) 
is temperature independent, H c i of Eq. (|56|l is only 
w eakly depe nds on T, whereas H * vanishes as T — > T c as 
i/l — T/T c . The last two are to be compared with the 
bulk Hci and H* , the bulk characteristics of the Bean- 
Livingston barrier, both of which are linear in 1 — T/T c 
near T r M 



V. DISCUSSION 

The behavior of narrow rings is closely related to that 
of long strips. The field given in Eq. I|55f) was derived pre- 
viously in a study of long strips of width W, where Hi 
was called the vortex exclusion fieldi^ According to this 
theory, when a superconducting strip is cooled through 
T c in an ambient magnetic field H a , vortices should be 
excluded from the strip when H a is less than Hi, be- 
cause the free-energy then has a global maximum at the 
strip center (at v m in our notation). On the other hand, 
when H a is greater than Hi, there is a local free-energy 
minimum at v m , where the probability of finding a vor- 
tex is proportional to exp(— V^ in /kBT). Near T c , this 
probability is close to unity. However, as T decreases, 
the characteristic energy scale, 0§/87r 2 A(T), increases 
rapidly. The vortices are then "frozen in" a potential 
well similar to those shown in Fig. [3] at a characteristic 
freeze- in temperature Tf, which for most superconduc- 
tors is estimated as being very close to T c . The value of 
A(Ty) turns out larger than strip widths of the order of a 
few /xm, which justifies calculations done in the thin-film 
limit. The above arguments suggest that SQUIDs to be 
cooled and operated in the earth's magnetic field should 
be made entirely of narrow lines in order to avoid flux 
noise due to thermal agitation of vortices trapped in the 
lines. Experiments by Dantsker et al. 22 are consistent 
with this conclusion. 

The solutions presented here for currents and energies 
in thin mesoscopic superconducting rings are of impor- 
tance for the physics of isolated rings (such as in Ref. 



.30, as well as for understanding the behavior of large en- 
sembles of interacting ringsi^ The statistical mechanics 
of these systems requires knowledge of ring energy levels 
and probabilities of quantum or thermally activated tran- 
sitions between the states of different energies. Our work 
should be useful for such statistical modelling. Despite 
the known shortcomings of the London approach, this is 
the only method that is of practical use for temperatures 
away from the critical temperature. A key advantage 
of the London equations is their linearity which makes 
exact solutions for mesoscopic rings possible. 

The electrostatic analogy we employed, may prove use- 
ful for various mesoscopic sample shapes. A number of 
these shapes (as squares, rectangles, or polygons) can be 
found in textbooks on applications of the theory of com- 
plex functions to the two-dimensional electrostatics, see, 
e.g., Ref. ^| An example of a thin-film disk is considered 
in the Appendix. 
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APPENDIX A: THIN FILM DISK 

It is worth stating upfront that the solution for a disk 
cannot be obtained by setting the inner radius a = in 
the solutions for a ring. The ring topology differs from 
that of a disk, and there is no continuous transition from 
the one to another. In fact, for a disk with no vortices 
the vorticity N = 0. The state of a disk is determined by 
the continuous variable H, the applied field, and by the 
vortex position (or positions, if any). 

Consider a thin-film disk of a radius b <C A with a 
vortex at an arbitrary position v. To find the current 
distribution one has to solve Eq. (JHJl with the boundary 
condition G(r = b) = const = Gh, or alternatively, Eqs. 
iJSJ for G v and l|10|l for Gh under boundary conditions 
G v (b) =0 and G H (b) =G b . 

The part Gh reads: 



G H (r) 



cH 
8^A 



r 

G In—, 

ro 



(Al) 



where the constants Go and ro are to be determined using 
the fact that the total current is I = —Gb- Evaluating 
/ with the help of London equation (|15|) one should take 
the phase change on circles r < v as zero, whereas for 
r>was — 2ir. Doing this algebra we obtain: 
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and 



G H {r) 



cH 
8^A 



In - 
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(A3) 



The problem of finding the part G v is equivalent to one 
for the electrostatic potential of a linear charge ccf>o /8ir 2 A 
at the point v parallel to the grounded metallic cylinder 
of a radius b. The solution is given, e.g., in Ref. Il9t 



where h = irb 2 H/(f>o (note that the field normalization 
here differs from the main text). In the vortex absence, 
J-q = (0q/327t 2 A)/i 2 , so that the barrier for the vortex 
entry is 
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G v = 



c4>o I b / r 2 + v 2 — 2rv cos ip \ 

ml - 4 / — - — I , ( A4) 
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v V r z + x z 



2rx cos ip 



It is now easy to find the field at which a local minimum 
first appears in the disk center (cPV^/dr 2 — at v = 0): 



where r, (p are cylindrical coordinates, the zero azimuth 
is taken as that of v, and x = b 2 /v. As discussed, the 
vortex self-energy e v = 4*o\G v ((p — 0, r — > v)\/2c is loga- 
rithmically divergent, so that we set r = v + £ and <y9 = 
to obtain: 
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The following is obtained in the same way as in the main 
text. The energy in zero field is 



E(0,v) = e v (v) + 
and the magnetic moment is 



8tt 2 A v 



A 4 = 



Mb 2 ~v 2 ) Hb 4 
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(A6) 



(A7) 



Hi = 



nb 2 



The field H cl is determined by F+(0) = 



rr <t>0 , 6 

7TO z 5 



(A10) 



(All) 



the result given by Fetter >i4 

Thus for H\ < H < H c ±, the vortex in the disk center 
is in metastable state, which becomes stable for H > 
H c i. Finally, the barrier to the vortex entry disappears 
altogether when the point vq at which — moves to 
the edge, v a — b — £: 



8tt 2 A 
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h 2 
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The thermodynamic potential T = E(0,v) — J^ 1 jidH 
follows: 



-T - (A8) 
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